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Q^' The effect of the penguin amplitude on extracting a from CP asymmetries in B' 







p'^p decays is studied using information on the SU(3)-related penguin amplitude 



D i in B~^ —J- K*^p'^. Conservative bounds on non-factorizable SU(3) breaking, small 

"^ ' amplitudes, and the strong phase difference between tree and penguin amplitudes, 

are shown to reduce the error in a in comparison with the one obtained using isospin 

K^ , symmetry in B — > pp. Current measurements imply a = [90 it 7 (exp)l'l5 (th)]°. 
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1. A major purpose for studying B and Bg decays is achieving great precision in 
Cabibbo-Kobayashi-Maskawa (CKM) parameters and providing precision tests for the 
Kobayashi-Maskawa mechanism of CP violation |Tj. CP asymmetries in properly chosen 
decays can be related with high precision to the angles (3, 7 and a = 7i — (3 — 'joi the 
CKM unitarity triangle |2lEllll. 

Currently the most precise single (hadronic-theory independent) determination of a or 
7 is based on the CP asymmetries Cl and 5*^ in B^ — > p'^p~ and on isospin symmetry in 
the B ^ pp system. The observation that the p mesons in B^ —>■ p^p~ are nearly entirely 
longitudinally polarized [SI El 13 IHj has simplified this study to becoming equivalent to 
an isospin analysis in B ^ inr [0] . Using isospin triangle inequalities pTH ITT] the current 
upper limit B{B -^ p°p°) < 1.1 x 10"^ |T^ implies a = (96 ± 13)° [3 EI, which includes 
an intrinsic error of 11° from the penguin amplitude alone, and only 7° originating in the 
measured CP asymmetries Cl and Sl. 

The purpose of this Letter is to suggest an alternative way for studying the penguin 
amplitude effect on measuring a in longitudinally polarized B^ -^ P^P- We relate the 
penguin amplitude in this process to the longitudinal amplitude in B^ -^ K*^p^ which 
is dominated by a AS* = 1 penguin contribution. The resulting error in a is shown to be 
smaller than in the isospin analysis of the three B —>■ pp decays. This will be argued to be 
the case in spite of a larger theoretical uncertainty caused by flavour SU(3) and further 
approximations entering the determination of the penguin amplitude in B^ -^ P^P~- 
The main point is that a large relative uncertainty in the penguin amplitude leads to 
only a small uncertainty in a, once the penguin amplitude is established to be small. 
Applications of flavour SU(3) to 5° -^ vr+Tr" and 5° -^ K+n- or B+ -^ K°n+ ^^ [H] 
involve a somewhat larger theoretical uncertainty in a because the ratio of penguin to 
tree amplitudes is considerably larger in B^ -^ tt+tt^ than in B^ -^ P^P- 

2. The amplitude for longitudinally polarized p mesons can generally be written as 

^i(5° ^ p+p-) = Te'^ + Pe'^ . (1) 

By convention T and P are positive, involving the magnitudes of the CKM factors V*fyud 
and V*fjVcd, and the strong phase 6 lies in the range — vr < 6 < tt. Time-dependence for 
longitudinal polarization is described in terms of two CP asymmetries Cl and Sl [IS| , 

r,(i?o(t)^p>-) + r,(i?o(t)^p+p-) = -^^ ^°^(^"^') + ^^ "^(^"^') • (2) 

The asymmetries Cl and Sl are given by 

1 - |AiP 2Im(Ai) 

/^ _ I PP' C _ ^ pp' fO\ 

^+\'^pp\ ^+\'^pp\ 

where 

Al{B^ - P+p-) ■ ^ ' 

Substituting dTJ into these definitions, one obtains 

_ 2r sin (5 sin(/? + a) 

1 — 2r cos 6 cos{/3 + a) + r^ 



sin 2a + 2r cos 5 sin(/5 — a) — r^ sin 2/3 
1 — 2r cos (5 cos(/5 + a) + r^ 

where 

r = ^ > (7) 

is the ratio of the penguin to the tree amphtude. 

In the absence of a penguin amphtude (r = 0) one has Cl = 0, Sl = sin 2a. For smaU 
values of r one finds 

Cl = 2r sin 6 sin(/3 + a) + 0{r^) , (8) 

Sl = sin 2a + 2r cos 6 sin(/? + a) cos 2a + O(r^) . (9) 



Given the value of /3 jHj, 

P = i2l.7tlir , (10) 

the two measurables Cl and Sl provide two equations for the weak phase a and for the 
two hadronic parameters r and 6. An additional constraint on r, S and a is needed in 
order to determine the weak phase. 

We will use the decay rate for a longitudinally polarized state in B^ -^ K*^p^ . The 
magnitude of the penguin amplitude dominating this process is related by flavour SU(3) 
to the magnitude of the penguin amplitude in B*^ -^ P^P~ ^E] • An additional constraint 
may, in principle, be obtained using the process B^ -^ K*^p~ for longitudinally polarized 
final states. In this case SU(3) relations apply to P, T and 6 and their SU(3) counterparts 
in B^ -^ K*^p^. However, so far only an upper limit has been measured for the decay 
rate of this process ^7] , and further information about the longitudinal fraction would be 
required. 

The amplitude squared for decays into longitudinally polarized K*^p'^ final states can 
be written as 




Al{B^ - i^*V)|^P_av. = FTT^ FP' = 21AFP' , (11) 



where fp = (209 ± 1) MeV and fx* = (218 ± 4) MeV are the vector meson decay con- 
stants JH]! and P is the penguin amplitude defined in (^. This equation defines a 
parameter F, which equals one when neglecting non-factorizable SU(3) breaking correc- 
tions (i.e. SU(3) breaking not in decay constants and form factors) in magnitudes of 
penguin amplitudes, and other contributions as discussed below. We now define a ratio 
of CP-averaged decay rates, 

^ ^ (MA.]' ^l{B^^K*^p^)+Tl{B-^K*^p-) 
\\Vcs\fK*) Tl{B^^p+p-) + Tl(b'~^P+P-) ' 

whose measurement provides a third constraint on r, 6 and a: 

-jl = . (13) 

1 — 2r cos 5 cos{(3 + a) + r'^ 



Eqs. (0), ® and (fTSj) give the three observables C^, Sl and TZ in terms of r, (5 and a. 
Assuming F is known permits a solution for a up to discrete ambiguities. 

3. We proceed to discuss the parameter F which, crudely speaking, relates the pen- 
guin amplitude squared in i? — > K*p to the one in i? ^ pp. At the amplitude level, 
the parameter F involves several effects. In addition to non-factorizable SU(3)-breaking 
it includes corrections from a color-suppressed electroweak penguin amplitude, penguin 
annihilation contributions [THl ITH] . and a doubly CKM-suppressed penguin amplitude. 
These corrections are usually thought to be small, so that F is expected to be near unity. 
We shall discuss each of the four corrections in turn. 

The neglect of non-factorizable SU(3)-breaking corrections is implicit in all applica- 
tions of SU(3) flavour symmetry to hadronic B decays. Within the current experimental 
uncertainties there is no evidence for the need of such a correction in the analysis of 
B decays to final states involving two pseudoscalar mesons {B -^ PP) [T31 CH d] and 
decays into a pair of pseudoscalar and vector mesons {B -^ VP) [20] • We assume that 
final states with two vector mesons {B -^ VV) are no different in this respect. In the 
QCD factorization approach [2I1I22]! non-factorizable SU(3) breaking corrections arise 
primarily from differences in light-cone distribution amplitudes of K* and p. This cor- 
rection is unlikely to exceed 15% at the amplitude level. The doubly CKM-suppressed 
penguin amplitude proportional to iKj'bKisl/lKfeV'csl ~ 0.02 is negligible, since no plausible 
dynamical mechanism is known which would enhance this amplitude without enhancing 
the dominant penguin amplitude. 

More important are the colour-suppressed electroweak penguin amplitudes in both 
B^ — i> p^p~ and B^ — >• K*^p'^, usually denoted Pew ^3 o'^ "^^ew CS], and a penguin- 
annihilation amplitude in B^ -^ P^P~, denoted PA jTE] or 2/3| ^H]- Since the dominant 
QCD penguin amplitude is smaller for B -^ VV than for B -^ PP, these two contri- 
butions are comparatively more significant in B ^ VV than in B ^ PP decays where 
they are often neglected. For orientation, a QCD factorization calculation of i? — > VV 
decays |2S] gives that the colour-suppressed electroweak penguin correction decreases F 
by about 0.1. The penguin-annihilation effect is about —0.3, and thus turns out to be 
the largest contributor to F — 1 in spite of being formally suppressed by l/rrib [22j. A 
global SU(3) fit to all B -^ VV decays, which requires more data, may eventually be able 
to check the size of penguin annihilation amplitude in these decays. One consequence 
of this contribution Ig,^ is a non-negligible branching ratio for longitudinally polarized 
Bs -^ P^P^ decays, on the order of a few times 10"''' |24j . 

A random scan through the input parameter space in the QCD factorization calcu- 
lation fHl I2n| that includes all four effects yields a nearly Gaussian distribution for F 
with F = 0.65 ±0.36. This estimate depends crucially on whether the annihilation model 
adopted in [221 predicts correctly the magnitude and sign of PA or 2^^. Since we would 
not like to rely on this assumption, we shall adopt the wider range, 

0.3<F<1.5. (14) 

Thus we are allowing a variation in P^ by a factor of five and in P by a factor larger 
than two. We will study below the sensitivity of the extracted error in a to this rather 
conservative range, showing that in spite of the large theoretical uncertainty allowed in 
F the determination of a is quite precise since data requires r to be small. 



4. We now describe the experimental status of Cl, Sl and IZ. The most recent 
measurements of Cl and Sl by the BABAR 7] and BELLE [Bj collaborations are 

_ f -0.03 ±0.18 ±0.09 , , . 

^ ~ 1 0.00 ±0.30 ±0.09 . ^ ' 



-0.33 ± 0.24+°-°^ 
0.08 ±0.41 ±0.09 



^ J -u.oo in u.^^_o.i4 , ,^„s 

^ 1 n na -I- n /I1 -I- n no v / 



Here (and below) BABAR and BELLE values are represented by upper and lower entries, 
respectively. These values imply the averages pKj . 

C^ = -0.03 ± 0.17 , ^L = -0.21 ±0.22 . (17) 

In order to compute TZ we use the CP-averaged branching ratios (given in units of 
10"^) and longitudinal polarization fractions, as obtained by BABAR jTj 123 I2Z] and 
BELLE [31211, 

E(n+n-^ _ / 30±4±5, , r 0.978 ± 0.014l°:°2i , 

^^^ ^ ^ - \ 22.8 ± 3.8+il , ^^^^ P )-'y o.941+0:01^ ± 0.030 , ^^^> 

,.r?<'*On+^ - /l'^-0± 2.9 ±2.01?:°, . .^*o,+^_/ 0.79 ±0.08 ±0.04 ±0.02, 
^^^ ^ ^ - j 8.9 ±1.7 ±1.2, ^^^^ ^ ^- j 0.43±0.1ll°:°t ^^^^ 

Using the B+ / B^ lifetime ratio t+/tq = 1.076 ± 0.008 25^, this implies 

r 0.0199 ±0.0065, , . 

I 0.0077 ±0.0032. ^^ 

The two values, representing BABAR and BELLE results, are not in good agreement 
with each other. The difference of 1.7cr originates mainly from a difference by a factor 
3.5 between the two measurements of longitudinal B^ -^ K*^p^ branching ratios. The 
weighted average of the two values in (PUJ) is 7^ = 0.0101 ± 0.0029. Calculating 7?. from 
the averages of (fTH)) and (fTU|) . we find a slightly larger value (implying a slightly larger 
error in the extracted value of a), 

7^ = 0.0125 ± 0.0031 . (21) 

We will use this value, the error of which does not include a scaling factor to account 
for the disagreement between the BABAR and BELLE measurements in (I19j) . We may 
expect this disagreement to disappear in the future. Note, however, that the effect of 
the experimental error in ()2H) on the extracted value of a is smaller than that of the 
theoretical uncertainty given by the wide range (|TH) for F to which 7^ is proportional 
[see (P|l ]. 

5. For given values of Cl, 5*^, TZ and fixed -^,(151),® and (fT!?|l can be solved 
numerically. The solutions exhibit an eightfold ambiguity for 6 and a in the range 



solution a[°] r 5 [°] (central) 

(1) 89.81^:? 0.124;0:°Jf -8 

(2) lOl.ete:? Clllt^M -172 

(3) n^.Stll O.lSltH^^ 28 

(4) 172.5;^°!^ 0.107lHf^ 152 

Table I: Four solutions for a, r and 5 corresponding to F = 0.9. 

— TT < 5 < TT, — TT < a < TT, wliicli cau bc understood and resolved into three independent 
invariance transformations (2'^ = 8) obeyed by (0), ® and p3p : 

(i) 5 ^ TT + (5, a -^ TT + a, r ^ r , 

(ii) (5^7r-(5, a^a[l + 0(r)], r-^r[l + C(r)] , (22) 

(iii) sin5^-sin(5^^^^^^i^[l + 0(r)], a ^ (37r/2 - a) [1 + C>(r)], 
cos(p — a) 

r ^ r [1 + 0(r)] . 

The first transformation is an exact symmetry of the three equations, leading to unphysical 
values of a larger than vr or negative. These four solutions can be discarded, leaving four 
solutions in the range < a < vr. The second and third transformations, 5 ^>- tt — 5 and 
a -^ 37r/2 — a, do not change the leading terms in r for C^, 5"^ and TZ given in (jH)), © 
and (fT^ . Including the non- leading terms in the expressions for C^, Sl and TZ implies a 
correction of order r in r, and corresponding corrections in a and 5, in the transformations 
(ii) and (iii). 

Keeping the theoretical parameter F fixed at its central value, F = 0.9, and using 
the measurements given in (fTTj) and (PT|) for Cl, Sl and TZ, we solve for r, 6, and a in 
the physical range < a < tt. The four solutions obtained within x^ = 1 contours for 
{Cl, Sl,TZ) are given in Table I. An important observation is the small value of r, in the 
range 0.10 — 0.13, which is implied by the small measured value of TZ. While the errors 
obtained for a and r are reasonably small, we only quote central values for S for which 
the errors are large. (See Figure 2 and discussion below.) We see that, as implied by the 
transformation (ii), the solutions (1,3) transform to the solutions (2,4) under 6 ^ n — 6. 
The change in a under this transformation is first order in r and is therefore rather small. 
A much larger change in a is implied by the transformation (iii), a — >■ 37r/2 — a, replacing 
(1) — *• (3), (2) — > (4). Solutions (3) and (4) are excluded by the measured value of P in 
(Uni) and by a + /3 + 7 = TT. 

The two remaining solutions, (1) and (2), both lying in the vicinity of a = 7r/2, 
can be distinguished by their values of the strong phase 6. It is clear from ©, where 
sm{/3 + a) cos 2a < holds for both solutions, that the smaller and larger solutions 
for a correspond to cos 5 > and cos 5 < 0, respectively. In the QCD factorization 
approach [211 122] the phase 6 is predicted to be small, being suppressed by I /nib or 
as{mb). This excludes solution (2) leaving as the single solution the value a = 89.8lll'l. 
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Figure 1: Dependence of a on F for the two solutions in Table I: (1) lower band, (2) 
upper band. The two bands denote experimental errors from Cl, Sl and TZ. 



Note that we do not need to assume that the phase 5 is small. It is sufficient to exclude 
values of 5 near ±180°. (A more precise requirement, depending on experimental errors 
on 5, will be given when discussing Figure 2 below.) The error ±7° in a is essentially the 
same as the error obtained in a^s using the isospin method [7]. This is not surprising, 
since by fixing the value of the parameter F to the central value in the range (fT^ we have 
restricted the origin of the error in the extracted value of a to experimental errors in the 
asymmetries Cl and Sl- The effect of the error in TZ given by (pT|l is relatively minor. 

The only theoretical error in our method (up to the discrete phase choice) originates 
in the parameter F. We now discuss the extraction of a for the entire range of F given 
in ()14|) . focusing on the two solutions (1) and (2) near a = 7t/2. In Figure ^ we show 
the dependence of these solutions on the parameter F. The lower and upper solid dark 
lines, corresponding to solutions (1) and (2) respectively, use central values for Cl, Sl 
and TZ. The bands around these two lines give experimental errors originating in these 
three measurements. Focusing on the theoretical error from F alone, we consider values 
of a along the dark solid lines, comparing values at F = 0.9 with values at F = 0.3 
and F = 1.5. We find the variation in the lower and upper solutions to be given by 
^-5!o)° ^^^ (101-6ti;2)°, respectively. We discard again the second solution on the 



basis of involving values of 6 in the neighborhood of vr rather than near zero. Including 
the experimental error from Table I and the above theoretical uncertainty from F, we 
conclude 

a = [89.8+^:? (exp)+^:^ (th)]° . (23) 

We note that the theoretical error, following from the range (|i^| ) in F and the preference 
for one of the two theoretically possible solutions, is considerably smaller than the error of 
11° in a obtained from an upper bound on \a—acs\ by applying the isospin triangle analysis 
to B -^ pp j7/. It is worth recapitulating the origin of this small error: Data on 7?. implies 
that the penguin correction is small. Once this is established the relation Sl = sin 2a 
receives only small corrections, and since sin 2a is rapidly varying near a = tt/2 even a 
significant error in Sl translates into a small error in a. 
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Figure 2: x^ = 1 contours for C^, Sl, IZ projected on the (a, 5) plane. Dotted, solid 
and dashed curves correspond to F = 0.3, 0.9, 1.5. The six points at focus denote the 
solutions for vanishing errors in C^, Sl, IZ. 



In order to quantify the criterion excluding solution (2) for a, we study now the 
dependence of the two solutions for a near 7r/2 on the strong phase 5. In Figure El we plot 
the x^ = 1 contours for Cl, Sl and TZ projected on the (a, 6) plane. Three different values 
oi F, F = 0.3, 0.9, 1.5, are described by the dotted, solid and dashed curves, respectively. 
The upper and lower parts of the curves correspond to solutions (1) and (2) discussed 
above. The x^ = 1 contours of the two solutions merge because of current experimental 
errors. The six points of focus for the three curves, marked on two almost parallel solid 
segments around 6 = —8° and 6 = —172°, are obtained for vanishing experimental errors 
in Cl, Sl and TZ. The length of the segments gives the purely theoretical uncertainty in 
a originating in the range of F. Figure El shows that when including current experimental 
errors in Cl, Sl and 7Z the two solutions (1) and (2) are presently distinguishable by 
\S\ < 7i/2 and \6\ > 7i/2, respectively. The additional requirement \6\ < tc/2, which 
excludes the second solution, will be relaxed considerably with more precise data on Cl, 
the error of which determines the uncertainty in 6. 

6. We conclude with a few comments about future improvements in the determination 
of a. 

1. The theoretical error in the extracted value of a depends weakly on the range as- 
sumed for the parameter F and on the measurement of TZ. A resolution of the 
disagreement between the BABAR and BELLE measurements of TZ will be reassur- 
ing. More precise measurements of Sl will have direct impact on the experimental 
error ±7° on a, while more precise measurements of Cl will eventually reduce the 
phase assumption to a discrete choice. This may be compared to the isospin-method 



for extracting a from B -^ pp, where further reduction of the error depends on what 
values the branching fractions will take. An intrinsic theoretical uncertainty at a 
level of a few degrees, caused in the isospin method by an / = 1 final state originat- 
ing from the p width [SHj and hy p — uj mixing [20] , niay potentially be resolved by 
studying with very large statistics the dependence of i? ^ 47r decay distributions 
on the invariant masses of pairs of pions near the p mass. 

2. Our suggestion for improving the determination of a replaces the application of 
isospin bounds in B ^ pp by theoretical input on the rough magnitude of F and 
a weak assumption about the the relative strong phase between the penguin and 
tree amplitudes in 5° — >■ p^p~. Currently the assumption \S\ < 7i/2 is required, but 
a weaker condition will suffice in the future. One possible test of this assumption 
consists of comparing globally the pattern of tree-penguin interference in B -^ pp 
and B -^ K*p decays. 

3. Information about a is also obtained from B^ -^ p^vr^ decays, which involve two 
ratios of penguin-to-tree amplitudes, in B^ —>■ p+vr^ and B^ -^ p^7r+. SU(3) ar- 
guments relating these decays to B -^ Kp and B -^ K*7t [31j, and a calculation 
based on QCD factorization ^H] show that these two ratios are small, in the range 
0.1 — 0.2, being on the smaller side in the second approach. The small ratios imply 
a small deviation of a from the value of aes obtained in the absence of penguin 
amplitudes ^|. Current data for time dependence in B^ — > p^vr^, given in terms 
of four observables, C,S,AC,AS |2S], imply a^s = (94 ± 4)°. An SU(3)-derived 
bound on the effect of penguin amplitudes, \a — a^fil < 9°, implies a = (94 ± 10)° 
when adding theoretical and experimental errors in quadrature PH |S2] . A more 
precise determination, a = (94 ±7)°, follows from the observable S alone using a 
QCD factorization calculation for amplitudes and strong phases P^F . Both determi- 
nations require stronger assumptions than those made in this work. However, the 
consistency of the most precise measurements of a (hence, 7) is impressive, allowing 
us to conclude that a is in the vicinity of 90° within a few degrees. 

Acknowledgements. This work was supported in part by the Israel Science Foundation 
under Grants No. 1052/04 and 378/05, and by the German-Israeli Foundation under 
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Note added in proof 

Shortly after this paper was submitted BABAR presented the new values B{K*^p^) = 
10.0 ± 1.7 ± 2.4 and fL{K*^p^) = 0.53 ± 0.10 ± 0.06 P], so that the BABAR and 
BELLE results are now in very good agreement [see p9|) ]. The new value of TZ equals 
0.0080 ± 0.0023 instead of (PT|) . This leads to the following changes in our results: Since 
7^ and F enter our analysis only in the combination TZ/F [see (J13j) ]. the new value of 7^ 
and -F = 0.9 is equivalent to the old value of TZ and F = 1.41 in Figure ^ It can be seen 
that within current experimental errors the two solutions corresponding to (1) and (2) in 
Table m overlap even for F = 0.9, i.e. the corresponding x^ = 1 contours in Figure |21 merge 
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also for F = 0.9. Separating the two solutions by requiring 6 < 7r/2, our final result (|^ 
now reads 

a = [91.2lli {exp)tli (th)]° . 

The larger experimental error is due to the fact that the two solutions have merged. 
As discussed in the text, the forseeable improved measurement of Cl will remedy this 
problem. 
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